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Abstract: We model-independently analyze the angular distributions of B — > Kll decays, 
I = e, fj,, for low dilepton mass using QCD factorization. Besides the decay rate, we study 
the forward-backward asymmetry A l FB and a further observable, FL, which gives rise to 
a flat term in the angular distribution. We find that in the Standard Model F l H oc mf, 
hence vanishing Ffj and F 1 ^ of around 2% (exact value depends on cuts) with a very small 
theoretical uncertainty of a few percent. We also give predictions for Rk, the ratio of 
B — * Kp,fi to B — > Kee decay rates. We analytically show using large recoil symmetry 
relations that in the Standard Model Rk equals one up to lepton mass corrections of the 
order 10 -4 including a s and subleading 1/E power corrections. The New Physics reach of 
the observables from the B — > Kll angular analysis is explored together with Rk and the 
B s — * 11 and B — > X s ll branching ratios for both I = e and I = fx. We find substantial 
room for signals from (pseudo-) scalar and tensor interactions beyond the Standard Model. 
Experimental investigations of the B — * Kp,fi angular distributions are suitable for the LHC 
environment and high luminosity B factories, where also studies of the electron modes are 
promising. 
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1. Introduction 



The exclusive decays B — > Kll with I = e,fi are governed in the Standard Model (SM) 
by flavor-changing neutral currents, and hence constitute sensitive probes of New Physics 
(NP). The three-body decays allow to study non-trivial observables by kinematical mea- 
surements of the decay products. They give access to a double differential decay spectrum 
with respect to the invariant mass of the lepton pair q 2 and a lepton charge asymmetry an- 
gle cos 6. In the absence of large statistics, partially integrated spectra such as the dilepton 
mass spectrum dTi/dq 2 or the angular distribution dTi/dcosO can be explored. Further 
= T(B — ► Kll) is in general different for electrons and muons. Having various theoretical 
or experimental advantages, the B — ► Kll observables cover a wide range of SM tests and 
NP searches, that are well suited for experimental study at high luminosity facilities at the 
T(4S) and the Large Hadron Collider (LHC), e.g., [1]. 

The B — > Kll branching ratio has been determined experimentally to be in agreement 
with the SM within uncertainties, and lies in the 10 -7 region [2, 3, 4, 5, 6]. Early data 
on more elaborate observables and q 2 -spectr& are beginning to come from the B factories 
[2, 3, 4]. While theoretical studies presented extensive phenomenological analyses of the 
dilepton mass distribution [6, 7], a detailed exploration of the SM background and NP 
potential of the angular dependence in the decay distribution is lacking. The B — > Kll 
angular distribution is very simple in the SM [6, 7] 

dT SM 

^— - ocsin^ + OK), (1.1) 
dcos9 

up to small lepton mass corrections of kinematical origin. A closer analysis shows that the 
cos 6>-dependence of the (normalized) angular distribution can be parametrized as [4, 6, 7] 



1 dT 



Tidcose =1(1-^(1 -cos 2 !?) + 1^ + ^0085, (1.2) 

with a flat term F l H /2 and a linear term in cos#, the forward-backward asymmetry ^4p B - 
Both are small within the SM, and therefore can signal the presence of NP. In particular 
they can be affected by Higgs and tensor interactions. Note that in the limit of vanishing 
lepton masses the SM predicts the same rates for electrons and muons r e SM = T^ SM + 
0(m£) if the same kinematical cuts are used [8]. 

In this paper we analyze the angular distributions of B — > Kll decays. We explicitly 
quantify the corrections to (1.1) within the SM and study model-independently the effects 
of (sb)(ll) operators induced by physics beyond the SM on (1.2). We use the framework of 
QCD factorization (QCDF) valid in the low-q 2 region [9, 10] and exploit the symmetries 
of QCD in the large recoil limit of heavy-to-light transitions [11, 12]. Also, resonance 
contributions from B — > K(cc) — > Kll can be controlled for dilepton masses below the 
charm threshold. 

The plan of the paper is as follows: After setting up the effective weak Hamiltonian 
in Section 2 hadronic matrix elements are given in Section 3. Section 4 contains model- 
independent formulae of the double differential and angular decay distributions. We give 
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numerical predictions for the SM in Section 5 including a detailed discussion of uncertain- 
ties. We also derive analytical expressions for T; and F l H obtained in the large recoil limit. 
In Section 6 we work out the sensitivity of the B — ► Kll angular distributions to NP in 
correlation with other observables in b — > sll decays. We summarize in Section 7. Techni- 
cal details about form factors and form factor symmetry relations in the \ow-q 2 region are 
given in Appendix A, whereas details on the B — > Kll hadronic matrix element in QCDF 
can be found in Appendix B. 



2. The Effective Hamiltonian 

The AB = 1 effective Hamiltonian 

n eS = -^vM^CiMOiM (2.i) 

is given in terms of dimension six operators Oi and their respective Wilson coefficients 
Cj. Both depend on the renormalization scale /j,, for which we take a low energy scale 
of the order of the 6-quark mass when evaluating i?-physics amplitudes. In (2.1) the 
leading CKM elements Vi m are factored out. The sum over i comprises the current-current 
operators i = 1,2, the QCD-penguin operators i = 3,4,5,6, the photon and gluon dipole 
operators i = 7, 8 and the semileptonic operators i = 9, 10. They are defined as 

2 
2 

° 8 = J^2 m b[s^ v P R T%)Gl v , 1O = j^[s^P L b][h^ll (2.2) 

where Pr/l = (1± 7s)/2 denote chiral projectors and m&(/x) the MS 6-quark mass at the 
scale /j,. For the operators Oi with i = 1, . . . , 6 we use the definitions given in [13], also used 
by [9, 10]. This set of operators suffices to describe b — > sll induced processes in the SM, 
which are dominated by Cj,C§ and Cio, whereas Cg enters at higher order in the strong 
coupling. 

Beyond the SM, NP might contribute in various ways. Assuming that NP manifests 
itself at and above the electroweak scale, it can be model-independently analyzed in the 
effective theory framework by allowing for NP contributions to the Wilson coefficients of 
the SM operators and by additional operators not present in the SM. To account also for 
the latter we include the most general b — > s (pseudo-) scalar and tensor operators with 
dileptons into our analysis: 



° lp = J^Y [ " m [h5l] ' = Sk? mb] [h5l] ' 



e 2 , .„ T ,.„., „, e 2 



&T = J^f [ * a ^ b] ' &T5 = [ * a ^ h] [Ia " Ul5l] ' (2 - 3) 
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where we made the dependence on the lepton flavor explicit by the superscript I. Note 
that there are only two independent tensor operators in four dimensions. At higher order 
also 4-quark operators with scalar, pseudoscalar and tensorial structure contribute to rare 
radiative and semileptonic decays [8, 14]. Here we neglect these effects. 

The additional NP operators (2.3) mix under QCD only with themselves. Their 1-loop 
anomalous dimensions 7 j = ^7^ are 

7 (°) = -6C F = -8, i = S,S',P,P', 

7 (°) = 2Cf = |, i = T,T5. (2.4) 

In our NP analyses all Wilson coefficients are taken at the low scale 



3. The Hadronic Matrix Element at Large Recoil 

A systematic treatment of the matrix element M. [B — > Kll] = {l(p^)l{p + )K{px)\H c s\B{pB)) 
is available in the large recoil region. We denote by Pb,Pk,P- and p + the 4- momenta of 
the S-meson, kaon, lepton / and antilepton I, respectively, and Mb, Mr- and mi are the 
corresponding masses. At large recoil the energy E of the i^-meson is large compared to 
the typical size of hadronic binding energies Aqcd *C E and the dilepton invariant mass 
squared q 2 = (p_ + p+) 2 is low, q 2 <C Afg. Consequently, in this region the virtual photon 
exchange between the hadronic part and the dilepton pair and hard gluon scattering can 
be treated in an expansion in 1/E using either QCDF or Soft Collinear Effective The- 
ory (SCET) [15]. Furthermore, only one soft form factor £p(q 2 ) appears in the B — > K 
heavy-to-light decay amplitude due to symmetry relations in the large energy limit of QCD 
[11, 12]. Other nonperturbative objects present are the light-cone distribution amplitudes 
(LCDAs) of the B- and i^-mesons, leading to numerically smaller contributions. This 
framework has been previously applied to B — > K*Tl decays using QCDF [9, 10] or SCET 
[16]. In this work we use the results from QCDF valid at low q 2 [9, 10] and include effects 
of finite lepton masses in B — > Kll decays. 

The B — > Kll matrix element can be written as 

M[B ^ Kll] = i^-V th V* s t P (q 2 ) {f vV % [l^l] + F A p% [l W l] (3.1) 

+ (F s + cos 9 Ft) [Tl] + (F P + cos 6F T5 ) [l l5 l]j . 

Here, 9 denotes the angle between the direction of motion of the B and the negatively 
charged lepton I in the dilepton center of mass frame, following [7]. Note that this conven- 
tion differs from other works, e.g., [6], where 9 is defined with respect to I. The functions 
Fi = Fi(q 2 ), i = S, P, A, V, T, T5 are given as 



2V%gz Mq 2 ) rl „ _2V\p i _ f T (q 2 ) . 

M B + M K f + {q 2 ) M B + M K f + (q 2 ^ 
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2 m b -m s f + (q 2 ) 



Mi 



Ml 



HQ 2 

Mq 2 



(3.2) 



F S 



IM% 



Ml 



fo(q 2 



2 m b -m s f+(q 2 ) 



(C l s + Cl), F V = C 9 + 



2m b T P (q 2 



+ 



8m/ 



Mq 2 ) rl 



M B £ P (g 2 ) M B + M K f + (q 2 ) 



where 



M% 



+ M i K + q i 



2(M' B M'l + M' B q 2 + M 2 K q 2 ^ 




-f, (3.3) 



and it is useful to note that 2p B ■ (p + - p_) = y/\p t cos 9. In the SM holds F| M = F| M = 
Fj^f = 0. Above, we have written the matrix element with the form factor £,p(q 2 ) = f+(q 2 ) 
as an overall factor. It constitutes the main source of theoretical uncertainties. The form 
factor ratios /o//+ and fr/f+ are constrained by symmetry relations at large recoil [11, 12], 
which are given in Appendix A together with definitions of the form factors and a discussion 
of their uncertainties. The quantity Tp(q 2 ) appearing in the vector coupling to leptons, 
Fy, takes into account virtual one-photon exchange between the hadrons and the lepton 
pair and hard scattering contributions. Tp(q 2 ) can be extracted from [9] and is given in 
Appendix B. At lowest order (denoted by the superscript (°)) up to numerically small 
annihilation contributions, it reads as 



4°\q 2 )=Zp(q 2 ) 



cru^xq 2 ) 



(3.4) 



Hence, Tp(q 2 ) takes care of the contributions from the 0\ q matrix elements oc Y(q 2 ) 
that are commonly included in an effective coefficient of the operator Og [17]. The next- 
to leading a s -corrections to Tp are known, see Appendix B, and taken into account in 
our analysis. Here we consider only NP effects from the NP operators (2.3), that is, their 
respective coefficients as appearing in (3.2) being non-zero, and Tp is SM-like. The 6-quark 
mass in Tp and Fy is the potential subtracted (PS) mass m^ s (/xj) at the factorization scale 
A*/ ~ \J ^QCD m i) and is denoted by m b throughout the paper. The 6-quark mass factors in 
F$ and Fp stem from the equations of motion, and we take them in the PS scheme as well. 
In the evaluation of the function Y(q 2 ) we use the pole mass mP h ° le [9]. The relation to the 
PS mass is given as m p b ° le = m^ s (/x/) + 4a s ^//(37r) [18]. The SM Wilson coefficients C 9 
and Cio are taken in NNLL approximation [13, 19]. The remaining SM Wilson coefficients 
Ci,...,6 and with their effective counterparts Cy S 8 enter only through Tp. For details see 
Appendix B and [9]. Note that chirality nipped operators O' 7910 can be readily included 
in the matrix element of B — ► Kll decays by replacing Cg^o — > Cg^o + Cg 10 in (3.2) and 
C 7 ->■ C 7 + C 7 in T P . 
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4. Decay Distributions of B — > 

Based on the matrix element (3.1) the double differential decay rate with respect to q 2 and 
cos 6 with lepton flavor I reads as 

''' ' ' = ai(q 2 ) + h(q 2 ) cos 6 + Cl {q 2 ) cos 2 9, (4.1) 



dq 2 cfcos 6 

where 

a ^(^ 2 ) _ 2 /«2| 77 |2 , I 77 |2\ , A 



g 2 (^2 + 1^2) + ^12 + ^1 



2\ 



+ 2nn{M 2 B -M 2 K + q 2 )Re(F P FX) + Am 2 M 2 B \F A \ 2 , (4.2) 

~ ^ , 2 = 2 {g 2 [PfRe(FsF^) + ite(F P i$ 5 )] 
1 o V A pi 4 p ^ 

+ mi [VXftReiFsF^ + (M 2 - M 2 K + g 2 )i?e(F T5 Fl)] }, (4.3) 
n(ql ' - q 2 (f3 2 \F T \ 2 + |F T5 | 2 ) - jP?(\F A \ 2 + \F V \ 2 ) + 2m l ^/xp l Re(F T F^) (4.4) 



and 

r ° " 512vr5M| • (45) 

These relations simplify considerably in the SM, where bf M (q 2 ) = and in the limit mi — > 
further holds af M (g 2 ) = -cf M (g 2 ). 

With (4.1) at hand the angular distribution 

d] 1 Ai + Bi cos 9 + Q cos 2 61 (4.6) 



cicos # 

is given in terms of the c/ 2 -integrated coefficients 



r,2 „„2 



flm&x f'max f max 

/ dg 2 aKg 2 ), 5/=/ ^ 2 6/(g 2 ), Ci=/ d? 2 Q(? 2 ). (4.7) 
J a 2 • J a 2 ■ J a 2 ■ 



'In 

Their values depend on the cuts in q 2 . We recall that while the boundaries of the phase 
space allow for dilepton masses in the rang e Am 2 <q 2 < (M B - M K ) 2 , our calculation is 
valid only in the low-g 2 region. Note that for very low dilepton masses there is sensitivity 
to light resonances. We therefore restrict our analysis to 1 GeV 2 < q 2 < 7 GeV 2 . 

The decay rate T/ and the integrated and normalized forward-backward asymmetry 
Ap B of the lepton pair can be expressed in terms of Ai,Bi and C\ 

T l = 2U l + 1 -Cy A 1 ¥B = y- ( 4 - 8 ) 
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We further introduce the observable 



2 . . „, r q 



2 



i i. .1 j ;■ i .n 



1 



Fi, = -(^+a) = y a d g 2 [ ai (? 2 ) + Q(? 2 )J / y 2 ^[^h^ 2 )]. (4.9) 

^min / ^min 

With (4.8) and (4.9), the angular distribution (4.6) is equivalent to (1.2) presented in 
the Introduction. Since F l H is normalized to T;, we expect reduced uncertainties in the 
former compared to the latter due to cancellations between numerator and denominator. 
As already anticipated after (4.5) within the SM a cancellation takes place in (4.9) between 
a/ and q such that Fjf M vanishes in the limit mi —> 0. We discuss this in detail in the 
next section. From here follows the approximate oc sin 2 9 angular dependence of B — > Kll 
decays in the SM as in (1.1). 

We would like to comment on the possibility of corrections to (1.2) or (4.6) from higher 
powers of cos 9, that is, a polynomial dependence in the angular distribution on cos n 9 with 
n > 2. Higher angular momenta arise from higher (> 6) dimensional operators in the 
weak Hamiltonian (2.1) or from QED corrections. Hence, they are suppressed by powers of 
external low energy momenta or masses over the electroweak scale, and a e /4ir, respectively. 
We discuss such corrections further at the end of Section 5 in the context of a non- vanishing 
forward-backward asymmetry in B — > Kll in the SM. 

A further useful observable in B — > Kll decays is Rk-, the ratio of B — > Kflfi to 
B — > Kee decay rates with the same q 2 cuts [8] 

= r ' F "- 4/3 °», (4.10) 

which probes lepton flavor dependent effects in and beyond the SM. We find that F l H and 
Rk are model-independently related 

R K -(1-F£-A) = 1, where A =\^^ - ^. (4.11) 

The expression for A simplifies in models where chiral couplings to electrons can be ne- 
glected as, for example, in the SM with m e = 0. Then Ffj = and T e = —4/3C e and in 
the SM A SM oc m 2 . We carefully examine SM predictions for T^F^ and Rk in Section 5 
and work out the NP potential of Rk, F l H and A l FB in Section 6. Corresponding values for 
A can be obtained by means of (4.11). 




5. Standard Model Predictions 

In this section we analyze B — > Kll decays within the SM. We give predictions for the 
observables F l H , Rk and Ti and the corresponding branching ratios B\ = B(B — > Kll) for 
low dilepton mass. Higher order SM contributions to the forward-backward asymmetry 
A l FB are briefly discussed. 

We start with a general analysis of lepton flavor dependence in the B — > Kll decay 
rate IY In the SM, such effects are of purely kinematical origin and often negligible [8]. At 
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large recoil, the suppression of the lepton mass induced terms can be quantified analytically 
using form factor symmetry relations (A. 3). For low q 2 , Ti can then be written as 



where the mj correction has been obtained from explicit expansion of the coefficients a\ 
(4.2) and q (4.4) in mj. (A useful relation is given in (A. 4).) Due to a cancellation with the 
kinematical function Pi there are no terms of order mf up to symmetry breaking corrections, 
which are estimated in the second correction term in (5.1). Form factor relations are broken 
in general by a s -corrections and power corrections in Aqcd /E, as discussed in more detail 
in Appendix A. As can be seen, these receive here further strong suppression from mf /Mg. 
Note that consistent with the Aqqd/E expansion we neglected terms of order M|-/Mj| 
and we approximated in the symmetry breaking correction contribution in (5.1) A « M|j, 
thereby dropping terms suppressed by q 2 /Mg. 

We conclude from (5.1) that lepton mass effects in the SM B — ► Kll decay rate at low 
q 2 are of order m^/q 4 ~ 10 -4 for muons and even further down by rn^/rn^ ~ 2 • 1CP 5 
for electrons, hence negligible in agreement with earlier numerical findings covering the 
whole dilepton mass region [8]. To leading order in mi, the decay rate depends then only 
on \Fy\ and \Fa\. The function Fa equals the Wilson coefficient Cf^ 1 ~ —4 with tiny 
dependence on the low scale Fy is a sum of Cg M ~ +4 and a term containing Tp. The 
latter is subject to unknown higher order power corrections. However, the typical order 
of magnitude of \Tp(q 2 )\ ~ 0.1 implies that these corrections constitute a rather small 
contribution to Fy and the corresponding uncertainties are very small compared to the 
dominating one from the overall form factor £p. 

The numerical analysis of T/ confirms the discussed qualitative features. The main 
uncertainties are due to the form factor £p, the CKM matrix element Vt s and the renor- 
malization scale ^b- For the form factor £,p{q 2 ) we use the findings from Light Cone Sum 
Rules (LCSR) [26]. At low dilepton mass, the form factor has an uncertainty between 
(12 — 16)%, with smaller uncertainty for larger q 2 , for details see Appendix A. Our nu- 
merical input is given in Table 1. We find that the /Uft-dependence of the decay rate is 
rather small, about a few percent, as can be seen from Figure 1 (left-hand plot). Here 
the coefficients a fl (q 2 ) and, to enable easier comparison, —c fl (q 2 ) are shown for Hb between 
m&/2 and 2mb. The small uncertainty due to \x\, is not unexpected because of the inclusion 
of NNLL corrections to the matrix elements of the current-current operators [27] in Tp, 
which cancels the ^-dependence of Cf M . In the right-hand plot of Figure 1 we show T M 
for three lower cuts q^ in = {0.5, 1, 2} GeV 2 as a function of the upper boundary g^ax- The 
combined uncertainty from £,p(q 2 ), lib and Vt s can be as large as 32%. Further subleading 
sources are the lifetime with 0.7% uncertainty and a e (/j,), which enters quadratically and 
brings in about 6% uncertainty to the B — ► Kll decay rates. The latter can be reduced 
by including the higher order electroweak corrections from [28, 29] to the renormalization 




(5.1) 
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a s (m z ) = 0.1176 ±0.0020 
a e (mb) = 1/133 
\V U \ = 0.0409 ±0.0021 
\Vd,\ = 0.0416 ± 0.0007 
m w = 80.403 GeV 



m 



pole 



(170.9 ± 1.8) GeV 



[20] 

[21] 
[21] 
[20] 
[22] 
[9] 



m b = (4.6 ±0.1) GeV 
m p c ole = (1.4 ± 0.2) GeV 
B(B -► XM) = (10.57 ± 0.15)% [20] 



fK = 

f By,,d 
fB s = 

af(l 



a 



a 4 

Ab,+( 

ep(o) 

r B ± = 
r B o = 
T B S = 



[20] 

[23] 
[25] 
[25] 
[26] 

1.5 GeV) = (0.458 ± 0.115) GeV [10, 24] 
= 0.327 ±0.053 [25,26] 
= (1.638 ±0.011) ps [20] 
(1.530 ±0.009) ps [20] 
(1.425 ± 0.041) ps [20] 



(159.8 ± 1.4 ±0.44) MeV 
= (200 ± 30) MeV 
: (240 ± 30) MeV 
GeV) = 0.06 ± 0.03 
GeV) = 0.25 ±0.15 
GeV) = -0.015 ±0.1 



Table 1: The numerical input used in our analysis. We denote by nib the PS mass at the factor- 
ization scale fif = 2 GeV. We neglect the strange quark mass throughout this work. 

group evolution which should capture the leading effect. For a complete higher order elec- 
troweak analysis the QED-corrections to the B — > K 11 matrix element should be calculated. 
In the corresponding calculation for inclusive B — > X s ll decays collinear logarithms of order 
a e /(4ir) ■ log (mj/m&) arise for low dilepton mass cuts [29]. The resulting splitting between 
electron and muon final states, however, diminishes after experimental cuts which separate 
electrons from energetic collinear photons. How much this matters for T e and T M and Rk 
cannot be answered until these QED-corrections are calculated. The uncertainties in Ti 
from the charm, bottom and top mass are 2%, 0.4% and 2%, respectively. 

In Table 2 predictions for the SM branching ratios of B~ — > K~p,fj, and B° — > K°flLi 
decays are given including the uncertainties from £p(q 2 ), Vf S and fib added in quadrature. 
The relative errors due to £p and ^ are given also separately. Lepton mass effects are 
negligible in Ti and B\ , and the decay rates and branching ratios with electrons agree within 
uncertainties with the corresponding ones with muons. The splitting of (9.5—9.7)% between 
the branching ratios of the B~ and B° mesons is dominated by the lifetime difference, but 
there is also a small isospin breaking contribution from spectator effects residing in Tp. 

In view of the insensitivity of T; SM to lepton mass effects for I = e, fx and with regard to 
its large form factor uncertainty it was proposed in [8] to investigate the ratio r^/T e , i.e., 
Rk (4.10). Our numerical analysis confirms a cancellation of the hadronic uncertainties in 
Rk also for low dilepton mass as can be seen in Figure 2. Here we show Rk for different 
cuts = {0.5, 1, 2} GeV 2 versus q^ illx - The combined uncertainty due to the form factor 
Cp(q 2 ) and the renormalization scale Lib is given by the bands and is tiny. This can be seen 
also from Table 2. The Lib and £p(q 2 ) induced uncertainties in Rk are of comparable size, 
of the order < 10~ 4 . The deviation of R^ from 1 is mainly due to the inclusion of effects 
of O {m^/q 4 } ~ 10~ 4 given in (5.1). Any measured deviation of Rk from 1 thus will signal 
NP which does not contribute equally to and r e as, for example, in the presence of 
non-universal lepton couplings. 

Similar to Rk also the angular observable F l H (4.9) is a ratio, where the overall factor 
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Figure 1: In the left-hand plot ai(q 2 ) and — ci(q 2 ) defined in (4.1) are shown for I = fi in the SM 
as a function of g 2 for the renormalization scale \ih between rrib/2 and 2mb- In the right-hand plot 
the SM B — > _ftT/2/x decay rate is given for three different cuts g min = {0.5, 1, 2} GeV 2 as a function 
of Q^ax- Here the bands take into account uncertainties from the form factor £p, /z& and Vt s . 



1.0006 




Figure 2: The ratio -Rr- in the SM for different cuts qf nin = {0.5, 1, 2} GeV 2 as a function of q. 
The uncertainties from the scale fib and the form factor are added in quadrature. 
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Tq (4.5) drops out and uncertainties can cancel. With the aid of the form factor symmetry 
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i.uuuou_ 00007 


+0.004 
-0.003 


+0.010 
-0.006 


1 00031 +0 0001 ° 

i.uuuoi_ .00007 


+0.004 
-0.003 


+0.010 
-0.006 


Rk 


i.uuu^r_ 00007 


+0.004 
-0.003 


+0.010 
-0.006 


1 nnn^8+ - 00011 

l.UUU^8_ .ooo07 


+0.004 
-0.003 


+0.010 
-0.006 


±.uuuoz_ 00007 


+0.004 
-0.003 


+0.010 
-0.006 


1 ooo33 +0 00011 

i.uuuoo_ .00007 


+0.004 
-0.003 


+0.010 
-0.006 




1 ooo39 +0 - 00011 

i.uuuoy_ 00007 


+0.004 
-0.003 


+0.010 
-0.006 


1 ooo4o +aooon 

l.UUU4U_o .00007 


+0.004 
-0.003 


+0.010 
-0.007 



Table 2: SM predictions for B^ (in units of 10 7 ), F£ and Rk for charged and neutral B-meson 
decays and different q 2 cuts (<^ in , q^ ax ) = (1,6), (2, 6), (1, 7), (2, 7) GcV 2 (from top to bottom). The 
uncertainties from the form factor £p(q 2 ) and the renormalization scale /ib varied between rrib/2 
and 2to& are also given separately in percent of the central value. The corresponding branching 
ratios with electrons, B e , agree within uncertainties with the ones with muons, B^. For details see 
text. 

relations (A. 3) and (A. 4) we obtain a simple expression for F l H in the SM at low q 2 : 

FjfM = 2m f JjL & &tf)Vj? m F A ? + \F V \ 2 ) (5.2) 

-L / J a 2 . Q 

where the denominator Ti SM is given in (5.1). From the lepton mass suppression of ai + q 
in the numerator of (4.9) follows F l ^ M oc mf, and i^ SM /F£ SM oc m 2 /m 2 such that Ffj SM 
is negligible. The cancellation between a^(q 2 ) and c^(q 2 ) is also visible from Figure 1 (left- 
hand plot). Note that the leading term of the integrand in the numerator of (5.2) is the same 
as the one in the denominator (5.1) except for an additional factor of (3i/q 2 — 1/q 2 - We 
therefore expect large cancellations of uncertainties in the ratio for low q 2 . This concerns 
the ones from the form factor, the renormalization scale, Vt s and unknown subleading l/E 
corrections in Tp. 

As expected the SM values of F^ are rather small, i.e., at the percent level, with 
the exact value depending on cuts. This can be seen from Figure 3, where F£ is shown 
for q 2 lin = {0.5, 1, 2} GeV 2 versus the upper integration boundary (? max . For g^ in = 
{1,2} GeV 2 , Fjj ranges between 0.015 — 0.05 depending on the values of g 2 nin and g max . 
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Figure 3: The observable Fg in the SM depending on qf nax for three cuts q^ in = {0.5, 1, 2} GeV 2 
(left-hand plot) and normalized to the central value (right-hand plot). The bands include combined 
uncertainties from and the form factor £,p(q 2 )- 

F£ becomes larger for smaller dilepton mass intervals and also for lower values of the lower 
cut <^ in . SM values of F£ are given for some low-q 2 cuts in Table 2. Within uncertainties, 
the predictions for B — ► K~ll and B° — > K°ll decays are the same. 

Indeed our numerical analysis of F£ exhibits strong cancellations of uncertainties. The 
form factor £p and m, induce uncertainties of comparable sizes of order one percent, see 
Table 2. The combined uncertainty from £,p(q 2 ) and //& is indicated by the small bands in 
Figure 3 and result in an < 2% uncertainty, see also Table 2. Power counting suggests an ad- 
ditional uncertainty from form factor symmetry breaking of order q 4 /Mg \J Aqcd JE ~ 3% 
in Ffl. We also allow for subleading power corrections to the hard scattering contributions 
at the order q 2 /M^a s ^ Aqqb/E ~ 3%, see Appendix A. Taking all this into account, 
F£ can be predicted with an accuracy of ~ 0(6%) in the SM, which is a high precision 
for an observable in exclusive .B-decays. Due to its huge suppression from m 2 , Ffj is a 
null test of the SM. Comparing our SM predictions for Rk and F^, the former is known 
even more precisely due to the cancellation of the 0(m^) -terms at leading order in a s and 
\/E and the stronger suppression of the symmetry relation breaking corrections in Ti (5.1) 
compared to the ones in (5.2). In order for the relation (4.11) to hold, A SM must be equal 
to _i?£ SM at the level of 10" 4 . 

As already mentioned in Section 4, operators in the effective theory of dimension 
higher than six or QED corrections induce additional contributions to the B — > Kll decay 
amplitude, which can modify the angular distributions. As for the higher dimensional 
operators, in the SM they are, for example, generated at one-loop by the Higgs penguin 
and the box with one charged pseudo Goldstone- and one VF-boson [30]. Contributions to 
scalar and pseudoscalar operators arise then at the order ~ mim^jm^. Plugging 
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this into (3.2), (4.3) and (4.8), a non-zero forward-backward asymmetry A 1 ^ 1 ~ mf/m^y 
is induced, which is too small to be experimentally probed. The corresponding SM tensor 
contributions have not been calculated, but they are subject to a similar 0{mimi ) /my V ) 
suppression, and negligible as the scalar ones in the B — > Kll observables. 

Higher order a e -corrections to exclusive B — > Kll decays have not been considered so 
far. Besides reducing the uncertainty from the overall a e (n) in the decay amplitude, radia- 
tive corrections can distort the decay distributions at the level of a e /(4-7r). The generation 
of an interesting cos ^-behavior from QED has been demonstrated for K — > nee decays. Ra- 
diative corrections via K — > ^77 enter the e^e-form factor in the matrix element, e.g., [31], 
which can be parametrized in our notation as Fy — > Fy + a e / (4tt) cos 9Fy, see (3.1). Note 
that Fy, Fy are functions of q 2 only. The extra power of cos 9 implies not only a non-zero 
A 1 ^ 1 ~ a e /(Ait)Fy /C$ M , but also a cos 3 0-term of order a e /(47r)iVcf M and a suppressed 
cos 4 #-term of order (a e /(47r)-FV) 2 hi the angular distributions. Unless the unknown cor- 
rection factor Fy is significantly enhanced (S> 1), it is unlikely that a e -corrections have 
observable consequences in B — > Kll decays. 



6. Beyond the Standard Model 

In the first part of this section we perform a model-independent analysis of B — ► Kll decays 
for I = e and / = /x. The size of the deviations from the SM in the observables Fjj, Rk and 
A l FB due to the NP operators (2.3) is estimated. We show this for four benchmark scenarios 
in Section 6.1 to Section 6.4. The second part of this section, Section 6.5, contains a brief 
discussion of models with (pseudo-) scalar and tensor operators. All NP Wilson coefficients 
are assumed to be real and are understood to be at the low scale /it, i.e., here C\ = C\{^). 
Leading logarithmic renormalization group evolution to the electroweak scale can be done 
with the anomalous dimensions given in (2.4). 

We start with some general considerations about the dependence of the B — > Kll 
observables on the NP Wilson coefficients. Up to corrections of order mf we find for 
the branching ratio 



1.64ps 



1.91 + 0.02 (Cf + Cj 2 ) + 0.06 (Cf + C l T \) + J£L - 



+ 



mt 



r<i2 

°T5 



C: 



12 



GeV 2 V3.28 2 3.28 2 10.36 2 5.98 2 



)+0(mf) 



(6.1) 



the numerator of F l H (4.9) 



2t b ± (Ai + Q) 



1.64ps 



mi 



( C T _ \ , m l ( 

V0.99 2.92/ GeV 2 V 



m. 



(0.51 GeV) 5 



+ 0.02 (C l s 2 + Cp 2 ) + 0.19 (C L 4 + C^) 



(6.2) 



C 



12 



r<i2 

°T5 



cj 2 



12 



GeV V 0.99 



GeV 2 V3.28 2 1.89 2 10.36 2 5.98 2 



) + O (mf) 



10" 
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and the numerator of the normalized forward-backward asymmetry (4.8) 



1.64ps 



om(c l s c l T + c l P c l Tb ) + 



GeW6.25 1.85, 



mi 



■( 



m f ( CgC? C l pC l T5 
GeV 2 V 4.12 2 4.12 2 



)+0{mf) 



10' 



(6.3) 



Here, we integrated over the dilepton mass region 1 GeV 2 < q 2 < 7 GeV 2 and used the 
central values of the input parameters given in Table 1. Then F l H is given by the ratio of 
(6.2) and (6.1), Rk by the ratio of (6.1) for I = /j, and / = e and A l FB by the ratio of (6.3) 
and (6.1), respectively. The contributions of the chirality flipped operators Og P can be 
included by the replacement C l s P — > C l s P + Cg P . 

As can be seen from (6.1), the B — ► Kll branching ratio is not very sensitive to NP 
effects from scalar and tensor operators due to the small coefficients in front of the NP 
couplings with respect to the large SM contribution. Moreover, the SM uncertainties of 
Bi will hide NP unless the Wilson coefficients become large, C' NP > 1. This actually can 
happen in some NP scenarios as we will show, in particular, in the decays into electrons, 
where the current experimental constraints are looser than the ones for the muons. Due 
to its tiny theory uncertainty the ratio Rk is a much more powerful probe of NP than the 
B — ► Kll branching ratios. Especially the terms at zeroth order in the lepton mass but 
also the ones linear in can significantly modify Rk — 1 with respect to its negligible 
SM value. 

The angular observables F l H (6.2) and A l FB (6.3) share several features with Rk — 1: 
They have a small and clean SM prediction and the sensitivity to tensor operators is 
higher than to scalar and pseudoscalar ones. Note that the dependence of B\ and F l H on 
the (pseudo-) scalar Wilson coefficients is the same and that the leading term in the lepton- 
mass expansion of A l FB requires the presence of both (pseudo-) scalar and tensor operators. 
Note also that Rk can be affected independently by NP in B —>■ Kee and B — > Kjlfj, decays. 

The available experimental information on F l H , Rk — 1 and ^4p B including SM predic- 
tions is given in Table 3 together with other related b — ► sll decay observables. The data on 
Rk include large dilepton masses where QCDF is not applicable and the ones on F l H and 
A l FB are in addition lepton flavor averaged. We do not take these constraints into account 
since they cannot be applied in a straightforward way besides having sizeable uncertainties. 

Important constraints on NP come from B(B S — > ll), which can be written as 



B{B S - ll) 



GWeM B J 2 B T Bs 

64vr 3 



\V tb V ts \ 



* |2 



Am 2 
Ml 



4mf 



m b + m s 



+ 



C P 



nib + m. 



C P 2mi 



Ml 



(6.4) 



The B s — > ll branching ratios depend on the difference of Wilson coefficients (C l s P — Cg P ). 
It follows that constraints from (6.4) can be evaded in the presence of both unprimed and 
primed (pseudo-) scalar Wilson coefficients unless there is a complementary constraint 
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observable 


sensitive to 


SM value 


data 




S,P °S,P> U T(5) 




0.8llg;g ± 0.46* [4] 


^FB 


S,P "+~ U S,P> U T(5) 


0(a e /(47r)) 


0.15^ ±0.08* [4] 
0.10 ± 0.14 ±0.01+ [3] 


R K - 1 


°S,P + U S,P> °T(5)' e VS - ^ 


C(10 J 


0.24 ± 0.31 1 [2, 4] 


B{B S /Z/x) 




(3.23 ±0.44) -10- 9 


< 8.0 • 10- 8 [5] 


- ee) 


^S.P ^S,P 


(7.56 ± 0.32) • 10~ 14 


< 5.4 • 10" 5 [33] 


[>0.04] 


U S IL, P ' U T(5) 


(4.15 ±0.70) • 10" 6 [6] 


(4.3 ± 1.2) • 10~ 6 [20] 


&T Cl [>0.04] 


°5 ^ °P ' L/ T(5) 


(4.15 ±0.70) • 10" 6 [6] 


(4.7 ± 1.3) • 10" 6 [20] 



Table 3: Observables in b — > sll induced transitions. Upper bounds are given at 90% C.L. For 
details see text. ^Data include g 2 -regions where QCDF does not apply and both I = e and \x are 
included. 



such as on (C l SP + Cg P ) from B — > Kll decays [8]. Tensor operators do not contribute 
to B s —>■ 11 decays and hence C l TT5 are not constrained by these decays. The current 90% 
C.L. upper bound on B(B S — > 11) for I = /x comes from CDF and D0 [5] 1 and for electrons 
from L3 [33]. The experimental information can be seen in Table 3 together with the SM 
predictions obtained with the input from Table 1. The bound on B(B S — > fifi) is O(20) 
away from the SM, and the one for electrons is nine orders of magnitude above the SM. 
As we show in Section 6.1, the current B(B S — ► ee) constraint is nevertheless on the verge 
of being useful, since NP in Cg P does not enter the B s — ► 11 modes with m/-suppression 
as the SM contribution, see (6.4). 

We further take into account the measurements of the branching ratios of the inclusive 
B -► A" s ee and B -► X 8 fifjt decays for q 2 > 0.04 GeV 2 denoted by #j ncl |[> .04]- The 
corresponding experimental values [20, 34] can be seen in Table 3 with SM predictions 
from [6]. The g 2 -cut dependent B — > X s ll branching ratios with (pseudo-) scalar and 
tensor interactions can be written as (see, e.g., [35]) 



B r\ [qLM - B(B -+ X s ll) = BrXlMsu + (\C^ + \C n ftMT (6.5) 
+ (\C l s + C l P \ 2 + \C y s + Cpf + \C l s - C l P \ 2 + \C y s - C l M 2 )M s , 



where 



M S,T = / <tf M 3 , T {q% B = 7' '^'g 1 ~7 ■ ': c ~ 1 ' (6.6) 

2m? J q i, (4tt)^ \V ch \ l f{m c /m b )K{m c /m b ) 



,2 



2 , 2^ « 3a 2 \V tb V*\ 2 B(B^XM 



mill 



and 



M s {q 2 ) = 2q 2 {m 2 - q 2 ) 2 , M T (q 2 ) = %im 2 - q 2 ?(2m 2 + q 2 ). (6.7) 



1 A stronger bound has been reported from a combined CDF and D0 analysis at 95% C.L., B{B S — » fifi) < 
5.8 • 1(T 8 [32]. 
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Here we neglect kinematical factors of m s and mi in the NP part and evaluate (6.6) and (6.7) 
with a 6-quark mass of 4.8 GeV, corresponding to the pole mass in accordance with [6]. The 
functions f(m c /mb) and K(m c /mb) represent the phase space function and QCD corrections 
of the decay B — ► X c li>i, respectively, and can be seen in [35]. Since M$,t > 0, NP from 
(pseudo-) scalar and tensor contributions enhances the B — ► X s ll branching ratios, and 
only the upper boundary of the experimental value of B(B — > X s ll) becomes a constraint 
on the corresponding Wilson coefficients. Also, since Mt 3> Ms, the inclusive branching 
ratios are more sensitive to tensor than scalar and pseudoscalar operators. Numerically, 
for 0.04 GeV 2 < q 2 < m\ we obtain M s = 1.92 • 10~ 8 and M T = 1.84 • 10~ 6 . 

In our NP analysis we also predict S™ cl |[ 1;6 ] for I = e,/x using Se nd |[i,6],SM = (1-64 ± 
0.11) • 10" 6 and B™\ lfi]jSM = (1.59 ± 0.11) • 10~ 6 [29]. These values are close to the 
experimental world average #i ncl |[i,6],ex P = (1.60 ± 0.51) • 10~ 6 [34] which is lepton flavor 
averaged and we therefore do not consider it as a constraint. However, we use this to 
illustrate the physics potential of future lepton flavor specific ^ e ncl |[i 5 6] and B™ cl |[i ; 6] mea- 
surements. The A^5,T-coemcients for this low dilepton mass region 1 GeV 2 < q 2 < 6 GeV 2 
are Ms = 0.52 • 10~ 8 and Mr = 0.83 • 10~ 6 . Note that we used here the 6-quark pole mass 
in the NP part of B(B -> X s ll) as well. To be consistent with the SM results of [29] the IS 
mass should be used once the next-to-leading order corrections to the NP part are known. 

Given the existing experimental constraints we cannot perform at present a fully model- 
independent analysis and fit for the six real NP Wilson coefficients per lepton species. 
Instead, we entertain in the following four benchmark scenarios with (pseudo-) scalar op- 
erators (Scenario I-III) and the tensor operators (Scenario IV) defined as: 

- Scenario I: NP in C l s and C l P , all other NP contributions vanish. 

- Scenario II: Same as Scenario I, but with the additional assumptions C l s = —C l p and 
C l oc mi. 

- Scenario III: NP in C l s , C l p and C§, C l p, the tensor coefficients C TT5 vanish. 

- Scenario IV: NP in the tensor coefficients C T , C T5 , all other NP contributions vanish. 

Scenario II is inspired by the Minimal Supersymmetric Standard Model (MSSM) for large 
values of tan (3, where tan/3 denotes the ratio of the two Higgs vacuum expectation values, 
see also Section 6.5. 

All B — > Kll observables are obtained for 1 GeV 2 < q 2 < 7 GeV 2 . We employ all 
bounds at 90% C.L. The resulting allowed ranges of the NP Wilson coefficients and the 
b — > sll decay observables are summarized in Table 4 and Table 5, respectively. Since the 
current experimental errors dominate the uncertainties, in the following we do not take 
into account SM uncertainties. Their inclusion would allow for slightly bigger NP effects. 

6.1 Scenario I: Scalars C l s and C l p 

We start with a discussion of the Wilson coefficients for muons, C$ p . The bounds on Cg P 
from (6.4) are displayed in the left hand plot of Figure 4, where contours of B(B S — ► fifi) < 
{0.05, 0.1, 0.2, 0.4, 0.6, 0.8, 1.0}- 10~ 7 are shown. The ranges for Cg p after applying the 90% 
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Wilson coefficient 


Sc I 


Sc II 


Sc III 


Sc IV 


fie 


[-8.3,8.3] 


— 


[-8.3,8.3] 


— 


C» 


[-0.69,0.69] 


[-0.55,0.41] 


[-5.6,5.6] 


— 


c% 


[-0.55, 0.82] 




[—5.6, 5.6] 




riel 






[-8.3,8.3] 




°S,P 






[-5.6,5.6] 




i~ie 








[-1.2,1.2] 


r<v 








[-1.1,1.1] 



Table 4: The allowed ranges for the NP Wilson coefficients C\ in Scenarios I-IV after using 
the constraints B(B S -> ee) < 5.4 • 1CT 5 , B(B S -» fifj,) < 0.8 • 1CT 7 , ^ ncl | >0 . 4 < 6.8 • 1(T 6 and 
BJ 1 ncl |>o.o4 < 6.3 • 10~ 6 , sec Tabic 3. A "— " means that the corresponding coefficient is zero in this 
NP scenario. 



C.L. B(B S — > Jifj) upper bound given in Table 3 can be seen in Table 4. The corresponding 
ranges of the observables are presented in Table 5. As can be seen, Fjj can deviate from the 
SM by about 40% whereas the forward-backward asymmetry is less then 1% in agreement 
with and updating earlier findings [7]. The deviation of the branching ratio B^ from the 
SM is less than 2% and completely negligible in view of the theoretical uncertainties. Also 
the NP contributions to 23™ cl |[ 1; 6] and £>™ cl |[>o.o4] are small compared to the theoretical 
uncertainties. 

The situation for the electrons is different due to the weaker bound from B(B S — ► ee) 
such that £>" lcl | [>o.o4] gives the strongest constraint on C s P . In the right-hand plot of Figure 
4 contours in the C% - Cp plane are shown for ^ ncl | [>0 .04] < {4.5,5.0,6.0,6.8,8.0} • 10~ 6 
and B(B S — > ee) < {0.1,0.5, 1.0} • 10~ 5 . The latter illustrates the constraints of improved 
measurements of B(B S — ► ee). We encounter the large ranges of Cg P given in Table 4 
allowed by £>" lcl |[>0.04] < 6.8 • 10~ 6 at 90% C.L., see Table 3. The corresponding ranges 
for the decay observables for I = e are presented in Table 5. As one can see, Ffj can be 
enhanced by orders of magnitude compared to its negligible SM value. Furthermore, Ffj is 
strongly correlated to B e , Rk (see (6.1) and (6.2)) and Hg ncI |r 16 i. The observables Ffj, B e 
and £>e ncl |[i,6] increase for increasing |C| P | whereas Rk decreases. We show B e versus Ffj 
(left-hand plot) and Rk versus Ffj (right-hand plot) in Figure 5. The branching ratio B e 
can be enhanced by about 60% with respect to its SM value. Bg ncl |[ li6 ] exhibits a similar 
enhancement but is subject to smaller theoretical uncertainties. The forward-backward 
asymmetry Ap B is negligibly small. 

The observable Rk depends on both lepton channels / = e, fi. In Scenario I the 
denominator B e receives large NP contributions whereas the numerator £> M stays close to 
its SM value due to the strong constraint from B s — ► fifi. This leads to a substantial 
decrease of Rk with respect to the SM as can be seen in Table 5 and also in the right-hand 
plot of Figure 5. 
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Figure 4: In the left-hand plot contours of B(B S — > /i/i) are shown in the Cg — Cp plane in Scenario 
I. The contours enclose values oiB{B s -> fifi) < {0.05,0.1,0.2,0.4,0.6,0.8, 1.0}- 10~ 7 starting with 
the innermost. In the right-hand plot contours of £>g ncl |[ >0 .04] < {4.5, 5.0, 6.0, 6.8, 8.0} ■ 10~ 6 (dashed 
black) and B(B S — ► ee) < {0.1,0.5, 1.0} • 10~ 5 (solid green) are shown in the Cf — Cf> plane in 
Scenario I starting with the innermost. 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 



F H F H 

Figure 5: Contours of Bf cl | [1)6] < {1.75, 2.0, 2.25, 2.35} ■ 10~ 6 in the F| - B e plane (left-hand plot) 
and the Ffj — Rk plane (right-hand plot) in Scenario I. 

6.2 Scenario II: MSSM-like C£ = -C£ 

Scenario II is a special case of Scenario I inspired by the MSSM in a certain limit (large 
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Observable 


Sc I 


Sc II 


Sc III 


Sc IV 


F e H 


< 0.39 


— 


< 0.56 


< 0.13 


F£ 


[0.013,0.035] 


[0.018,0.032] 


[0.013,0.56] 


[0.014,0.18] 


Rk 


[0.61,1.01] 


[0.996, 1.01] 


[0.44,2.21] 


[0.93,1.10] 


B e [10- 7 ] 


[1.91,3.14] 


— 


[1.91,4.36] 


[1.91,2.00] 


B, [10- 7 ] 


[1.90, 1.94] 


[1.90, 1.93] 


[1.90,4.26] 


[1.87,2.10] 


A% B [%] 


[-0.02,0.02] 


— 


[-0.02,0.02] 


[-0.02,0.02] 




[-0.6,0.6] 


[-0.5,0.3] 


[-4.46,4.46] 


[-3.1,3.1] 


B{B S - ee) [10^ 5 ] 


< 1.17 


— 


< 2.33 


— 


£(£ s -► ////) [10 '] 


< 0.8 


< 0.8 


< 0.8 




4 ncl l[i, 6] [io- 6 ] 


[1.64,2.35] 




[1.64,2.35] 


[1.64,2.83] 


4 ncl ki,6] [io- 6 ] 


[1.59,1.60] 


[1.59, 1.60] 


[1.59,2.17] 


[1.59,2.56] 


^ nCl |[>0.04][10" 6 ] 


[4.15,6.8] 




[4.15,6.8] 


[4.15,6.8] 


^ nCl |[>0.04][10" 6 ] 


[4.15,4.18] 


[4.15,4.17] 


[4.15,6.3] 


[4.15,6.3] 



Table 5: Allowed ranges for b — > sll observables in Scenarios I-IV after taking into account the 
constraints from B(B S — ► 11) and £>J ncl |[>o.o4] for Z = e and I = [i, see Table 3 and the text for details. 
A "— " means that the corresponding observable is SM-like. 



tan/3), see also Section 6.5. In this model, the (pseudo-) scalar Wilson coefficients are 
proportional to the lepton mass C l s p ~ mi, such that Cg P can be neglected and b — > see 
decays are SM-like. Furthermore the relation Cg = —Cp holds and the primed coefficients 
Cg p are suppressed by m s /m(, and can be neglected. 

The allowed range of Cg and the effects of NP on the rare decay observables are given 
in Table 4 and Table 5, respectively. Since Scenario II is a constrained variant of Scenario 
I the deviations from the SM are smaller in the former. The NP contributions to do 
not exceed 30% whereas the deviations of from the SM are of the order of 2%, much 
smaller than the theoretical uncertainties. The same holds for Z?™ 01 ^^], which confirms 
earlier studies within the MSSM [36]. Since B e is SM-like in Scenario II, the deviation of 
Rk from the SM is much reduced with respect to the one in Scenario I. We find NP effects 
of 1%, which are larger than the uncertainties of the SM prediction. The forward-backward 
asymmetry is smaller then 1% in agreement with previous works in the framework of the 
MSSM [37]. 

6.3 Scenario III: Scalars C l s , C l P and Cg, C l /> 

In Scenario III we use the full set of (pseudo-) scalar Wilson coefficients including the 
chirality flipped ones Cg P for I = e and / = fi. The constraint from the B s — > 11 branching 
ratios alone can be evaded due to cancellations between C l s p and C$ p , see (6.4). To 

obtain constraints on Cg P we combine B{B S — > 11) with jf5; ncl | [>o.04] data. We find the 
allowed ranges for the Wilson coefficients given in Table 4. In the electron sector C| P can 
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Figure 6: Contours of B" cl |[i l6 ] < {1.75,2.0,2.17} ■ 1CT 6 in the F£ - R K plane in Scenario III 
(left-hand plot). In the right-hand plot contours of Bf c \ xfi] < {1.75,2.0,2.25,2.35} • 10~ 6 are 
shown in the Ffj — Rk plane in Scenario III. For details see text. 

be as big as in Scenario I with identical ranges for Cap. In the muon sector the Wilson 

coefficients Cgp are now comparable in magnitude to the ones for electrons. 

The large Wilson coefficients lead to big NP effects in the rare decay observables, see 
Table 5. In Scenario III Rk can both increase and decrease significantly with respect to the 
SM as opposed to Scenario I where B(B S — ► /2/i) permits only a large decrease of Rk- The 
substantial deviations of Rk from the SM are already challenged by existing data given in 
Table 3. However, as already stressed, since these data contain also large-q 2 events where 
QCDF is not applicable it is not clear how to impose these constraints in a well-defined 
way. 

The increase of both F£ and Rk for increasing values of B]f cl | [i^] can be seen in the left- 
hand plot of Figure 6, where contours of £>™ cl |[ 16 ] < {1.75, 2.0, 2.17} • 10~ 6 are shown. Sim- 
ilarly, the increase of Ffj and decrease of Rk for increasing values of jSg ncl |[ 1)6 ] is displayed 
in the right-hand plot of Figure 6 with contours of ^ ncl | [1>6] < {1.75, 2.0, 2.25, 2.35} • 10~ 6 . 
The NP contributions enhance both B e and /?« by order 200% above the SM such that 
measurements of these observables in the low-g 2 region could provide constraints regardless 
of the large form factor uncertainties. Scenario III allows for |^4p B | < (4 — 5)% whereas 
^4| B is negligibly small. 

6.4 Scenario IV: Tensors C l T , C l T5 

In Scenario IV we consider only NP in C l T T5 . These Wilson coefficients do not contribute 
to B s — ► 11 decays and hence are currently constrained only by inclusive B — ► X s ll decays 
(6.5). The corresponding bounds can be seen in the left-hand plot of Figure 7, where 
contours of £>e nd |[>0.04] < {4.5,5.5,6.3,6.8,8.0} • 10~ 6 are shown in the Cf- — C|> 5 plane 




Figure 7: In the left-hand plot contours of £?g ncl |[ >0 4] are shown in the Cf. — Cf> 5 plane in 
Scenario IV. Each contour encloses values of #e ncl |[>o.o4] < {4.5,5.5,6.3,6.8,8.0}- 10~ 6 starting 
with the innermost. Corresponding constraints for — C^ 5 can be read off from the left-hand 
plot as well. In the right-hand plot contours of Sjf cl |[i, 6 ] < {1.75,2.0,2.25,2.56}- 10~ 6 are shown 
for Fjg versus Rk in Scenario IV starting with the innermost. 



starting with the innermost. The constraints on C^ T5 from upper bounds on £>™ cl |[>o.04] 
can be read off from the same plot. We find the ranges of the Wilson coefficients given in 
Table 4 using the 90% C.L. constraints B™ cl \ [>0 . 04] < 6.8 • 10" 6 and Bf cl \ [>0 . 04] < 6.3 • 1CT 6 , 
see Table 3. As anticipated after (6.7), i3| ncl |[ >0 .04] constrains C^ T5 stronger than Cg'p. 

The NP effects in and F^ are huge with respect to the SM predictions as can be 
seen in Table 5. F l H increases for increasing 0| ncl |n 6 i. This correlation is shown in the 
right-hand plot of Figure 7 for contours of £jf cl | [lj6] < {1.75,2.0,2.25,2.56} • 10~ 6 in the 
F£ — Rk plane. Similar correlations hold for the electron sector. Rk receives order 10% 
corrections from NP which are well above the theoretical uncertainties. The branching 
ratios Bi are subject to NP contributions < +10%, which cannot be separated from the 
larger form factor induced uncertainties. On the other hand, the NP enhancements due 
to C l TT ^ in S| ncl |[ 16 ] are larger, about 70%, which makes the inclusive decays a sensitive 
probe of tensor operators. As in all other Scenarios I-III is negligibly small. |^4p B | 
does not exceed 3%. 

6.5 Models with Scalar and Tensor Interactions 

While there are several known models beyond the SM with large (pseudo-) scalar inter- 
actions, tensor operators are often neglected. Let us begin with some general remarks on 
the origin of b — > sll tensor operators in the SM and the MSSM: In the SM they arise only 
at higher order in the electroweak operator product expansion (OPE) from finite exter- 
nal momenta in the matching calculation. In the MSSM tensor operators are induced at 
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leading order OPE only from photino and zino box diagrams, which are, however, sublead- 
ing in tan (3 with respect to the Higgs penguins discussed below. Higgsino contributions 
to tensors are further suppressed by down-type quark and lepton Yukawa couplings [38]. 
In addition tensors with two leptons are induced by scalar operators under QED renor- 
malization group running, hence are of higher order in a e /(47r) • log(/////&), e.g., [8, 14]. 
Another mechanism to generate tensor contributions is to consider models with scalars 
having appropriate quantum numbers such that tree level exchange induces the operators 
(jLbn)(sLlji) or (/r6i,)(sr/£,). Subsequent fierzing then leads to tensor operators. Among 
this class of models are those with leptoquarks. We consider such models below after briefly 
commenting on the MSSM at large tan/3 and the MSSM with broken i?-parity. 

For large values of tan (5 the MSSM produces substantial scalar couplings C l s P from 
Higgs penguins, for example, induced by chargino loops 

C^oc^tan 3 /?, 4',p-^4,p, (6-8) 

for exact Wilson coefficients see [7]. Here, m^o denotes the mass of the pseudoscalar Higgs 
boson. The relation C l s = —C l P holds only at leading order in tan/3 [7] and prevents the 
generation of tensor couplings from QED running [8]. The flipped coefficients Cg P are 
suppressed by the mass of the strange quark. Since C l oc mi the couplings to electrons are 
negligible. 

In i?-parity violating supersymmetry scalar and pseudoscalar FCNC-operators can be 
generated at tree level from the superpotential (see, e.g., [39]) 

W R = \ ijk V L L{e R + X^Li&A (6.9) 

where Ll(Ql) an d eR(dR) denote the superfields containing the lepton (quark) doublet 
and the charged lepton (down-type quark) singlet, respectively. One gets at the matching 
scale from sneutrino exchange 

r l „ ( 4? 2 A fc23 A fc» r U „ ( 4?r ) 2 X 'k32 X tll (ct 1 n ^ 

Cs ' P *—^G^f; ^^—VMGRmi- (fU0) 

Here denotes the sneutrino mass and summation over the sneutrino flavor k is under- 
stood. Contributions from squark exchange modify only the vector and axial vector type 
operators and are not shown. The couplings in (6.10) obey C l s = —C l p and C l g = +C l p . 
As in the MSSM at large tan (3 with unbroken i?-parity discussed previously, there are no 
tensor operators generated from leading order matching. The chirality flipped contribu- 
tions can be sizeable and help to escape the constraint from B s — > 11. The couplings (6.10) 
are then essentially only constrained by the B — > (K^*\X s )ee and B — ► (K^*\X s )p,fi 
branching ratios. The size of the possible modification of Rk from one is then given by 
the theoretical and experimental uncertainties of the branching ratios, see also [40]. 

The Lagrangian of leptoquarks <$lq coupling to a lepton and a quark can be written 

as 

£>LQ-l-q = \jkQj<t>LQ, (6-11) 
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where i,j label the lepton and quark generation, respectively. <$>lq can be a scalar or a 
vector under space-time transformations and a singlet, doublet or triplet under SU(2). 
Details can be seen, e.g., in [41], where also contributions to flavor-changing scalar and 
pseudoscalar operators have been discussed. Here we consider only the contributions to 
tensor operators, which are induced by tree level scalar leptoquark exchange and fierzing as 
explained earlier. The SU (2)-properties of the requisite operators require mixing of lepto- 
quarks with different »Si7(2)-quantum numbers. The latter is induced by interactions with 
the Higgs boson and arise after electroweak symmetry breaking, see the second reference 
in [41]. As a result, tensor operators in leptoquark models are suppressed by the vacuum 
expectation value of the Higgs v =< H° >= (2v / 2Gi?)~ 1 / 2 over the scalar leptoquark mass 
wig. Specifically for b — > sll transitions this yields the tensor coefficients 



Ct ' T5 " e 2 V th V t * s G F m% ml ' (fU2) 



7. Summary 



We thoroughly investigated the angular distributions in B — ► Kll decays in a model- 
independent way. We find that the cos ^-dependence in the normalized 1/T; dTi /cfcos 6 
spectrum, see (1.2), offers great opportunities to test the SM and search for NP. The req- 
uisite observables are the flat term in the distribution, F l H /2 and the forward-backward 
asymmetry ^4p B - The coefficient of cos 2 is related to F l H . No powers of cos greater than 
two appear in the B — > Kll angular distribution up to higher dimensional operators not 
present in TC e g (2.1)-(2.3) and QED corrections. Both are strongly suppressed by powers 
of the low energy masses and momenta over the scale of electroweak NP and by a e /(Air), 
respectively. 

In the SM, F l H oc mf, and Ffj is negligible. The SM value for F£ is small, order 
few percent, and can be cleanly predicted using QCDF for low dilepton masses with 2% 
accuracy, see Table 2. Taking into account subleading 1/S-corrections the uncertainty is 
conservatively inflated to ~ 6%. The forward-backward asymmetry vanishes exactly in the 
SM up to the aforementioned higher order OPE and QED corrections. The a e -corrections 
induce the parametrically leading contribution of the order a e /(47r). 

We also give SM predictions for the B~ — ► K'jlfi and B° — ► K p,jjL branching ratios. 
They have a substantial uncertainty of order 32% mostly from the form factor. On the 
other hand, the SM ratio of B — > Kjx\i to B — > Kee decay rates, R^, equals one at the 
level of 1CT 4 . We show analytically at large recoil using form factor symmetry relations that 
the apparent huge suppression of lepton flavor effects in results from the cancellation 
of C(m 2 )-corrections to leading order in 1/E and a s in the B — ► Kll decay rate, see (5.1). 
In addition potentially large corrections to Rk can arise from collinear QED logarithms, 
whose actual net effect depends on experimental cuts [29]. The corresponding calculation 
for B — > Kll decays has not been done. 

Beyond the SM, the observables F l H , A l FB and Rk are sensitive to Higgs and tensor 
interactions. We work out NP signatures and correlations by taking into account existing 
data on B(B S — ► 11) and B(B — > X s ll) for I = e and I = fi separately. We find that the 
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NP modifications to the angular observables Ffj, F^, Ap B and Rk — 1 can be sizeable, 
see Table 5. Even larger effects in the forward-backward asymmetries Ap B and Ap B arise 
in models where both (pseudo-) scalar and tensor operators are present. From a scan of 
twelve real NP coefficients Cg'p,C! rT5 for I = e and I = pL we find model-independently 
the upper bounds 

\A% B \ < 13%, \A£ B \ < 15%. (7.1) 

Both Rk and F£ enable precision tests of the SM in exclusive B — > Kll decays, but 
their experimental requirements are different: Whereas Rk requires only measurements of 
decay rates into both electrons and muons, Fj^ is extracted from the muon channel alone, 
however, at the price of an angular analysis. The latter needs high statistics and is well 
suited for the LHC(b) setup. NP searches with angular distributions in B — > Kfl/j, should 
also be feasible at the Tevatron, where CDF has recently measured B{B + — > K + p,fi) and 
B(B° -► K*°fin) [5]. 

The experimental situation for the observables F l H , A l FB and Rk is currently at a very 
early stage, see Table 3. In particular, all measurements average / = e and I = [i final 
states except the ones of Rk [2, 4]. Ultimately all observations in rare semileptonic decays 
B — > Kll,B — > if*ZZ and 2? — ► X S ZZ should be available for each lepton flavor separately 
since deviations from the SM could be /-dependent. For example, NP in the electron channel 
could escape the B — ► Kp,[i decay studies completely implying also Rk < 1- Existing data 
on the I = e modes are weaker than the corresponding ones for decays into muons, allowing 
for larger NP effects in the electron modes. In this way, b — > see induced channels such as 
B — ► ifee provide unique opportunities for the clean B factory environment. Appropriate 
cuts in q 2 should be taken into account to maximally exploit the theoretical predictions. 
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A. B — > Kll Form factors 

In this appendix we give definitions and properties of the heavy-to-light form factors for the 
B — > K transition at large recoil. The symmetry relations emerging in this region between 
the QCD form factors are reviewed including symmetry breaking corrections. Furthermore, 
details about the form factor /+ from Light Cone Sum Rules calculations [26] can be found 
here. 
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The B — > K matrix elements are parametrized in terms of the three QCD form factors 
f+Jo and fr as [42] 

q 2 



(K{p K )\- S1 ^b\B{p B )) = (2p B - q),U(q 2 ) + B g K q„[f (q 2 ) - f + (q 2 )}, (A.l) 



(K(p K )\sia^b\B(p B )) = -[(2p B - q),q 2 - (M 2 - 4)^ M ^ ' ^ 

At leading order in the 1/E expansion / +! o,t(<Z 2 ) obey symmetry relations [11, 12] such 
that they all can be related to a single form factor denoted by £,p(q 2 )- Within QCDF 
a factorization scheme has been chosen with f + (q 2 ) = £,p(q 2 ) [12]. Including subleading 
corrections, the symmetry relations can be written as 

/ n 2 

l + 0{a s ) + 



f+ M B 
fr M B + M K 



A 



QCD 



M% V E 



\ + 0{a s ) + 




(A.3) 



/+ M B 

up to higher order QCD, power and mixed corrections. The a s -corrections from the soft- 
overlap and hard scattering contributions indicated in (A.3) have been calculated in QCDF 
and are given in [12]. These corrections are taken into account in the numerical analysis 
of this work. Analogous relations can be found in the framework of SCET using MS 
subtractions [15, 43]. The symmetry relation breaking corrections due to subleading orders 
in the Aqcd/F expansion have been considered for the soft-overlap part using SCET [44]. 
The corresponding corrections are indicated in (A.3). Note that the expansion parameter 
is rather \/Aqcd/^ than Aqcd/F, and that for /o//+ an additional suppression of q 2 /M B 
appears. Subleading contributions from hard spectator scattering to (A.3) are unknown 
and arise at higher order, 0(a s ^Aqqd/E). 

The form factor symmetry relations (A.3) imply in the SM for B 

^\F P \ 2 + A\F A \ 2 + M2B -^ K + q2 2Re(F P FX) 

which enters af M + cf M , see Section 4. Here, the explicit SM expressions for Fy,A,p (3.2) 
have been used and F P = miFp has been rescaled. The relation (A. 4) involves only the 
ratio /o//+ and results in a beneficial q 2 /M B suppression of the power corrections. Beyond 
the SM the corresponding expression depends on all functions Fj and there are in general 
no cancellations from symmetry relations. 

We employ the form factor f+(q 2 ) = £,p(q 2 ) from LCSR calculations [26]. It is given 
in terms of the Gegenbauer moments of the K- meson LCDA, af and af as 

U(q 2 ) = f?(q 2 ) + a?( Pm )f?(q 2 ) + 4 (»m)ft 2 (q 2 ) + 4 (pm)fT(q 2 )- (A.5) 

The c/ 2 -dependent functions , i = 1, 2, 4 and f+ s are obtained from a fit and parametrized 
in [26]. Here we use "set 2" with m?°' 6 = 4.8 GeV corresponding to the infrared factor- 




ization scale fijp = \J M B — m^ ole2 = 2.2 GeV. The running of the Gegenbauer mo- 
ments given in Table 1 from 1 GeV to 2.2 GeV is accounted for by the scaling factors 
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Figure 8: The form factor £p(q 2 ) = f+(q 2 ) in the low-q 2 region including uncertainties from the 
Gegenbauer moments af (lighter shaded area) and from af and A as with their uncertainties added 
in quadrature (darker shaded area), for details see text. In the left-hand plot is shown £p(g 2 ), and 
in the right-hand plot the form factor normalized to its central value, £p(g 2 )/£p ntral (a 2 ). 

{0.793, 0.696, 0.590} for {af , af , af}. The relative uncertainty of /+ due to the asymptotic 
form factor f+ (which is independent of the af) at q 2 = is approximately A as // + (0) = 
10%, see Table 2 of [26]. In order to estimate the form factor uncertainty in the low-q 2 
region we scan over the Gegenbauer moments according to the ranges in Table 1 trans- 
lated to fijR = 2.2 GeV and add the uncertainty from A as in quadrature. The form factor 
f+(q 2 ) = £,p(q 2 ) with its uncertainties with and without A as is shown in Figure 8. The 
total uncertainty is 16% at maximal recoil and reduces to 12% at q 2 = 7 GeV 2 . The re- 
duction of the relative form factor uncertainty towards larger values of q 2 stems from the 
increase of the form factor in this region while keeping A as from q 2 = 0. The decrease of 
the form factor uncertainty for q 2 > has been considered likely in [26]. 



B. T P for B -> K 

The amplitude Tp{q 2 ) can be extracted from [9] as 



Tp = i P 

7T 2 f B fK 

N c M B 



C (0) + a s°F 
P 47T 



4 f) + 4 nf) 



(o) a s C F 



,(f) 



L P,± 



(B.l) 



where all Tp± , Tp ^ and T p n 2 are functions of (oj, u). fp and fx denote the B- and K- 
meson decay constants, respectively, whereas $>b,±{w) and are the corresponding 
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LCDA's. The remaining quantities are calculable perturbatively 

r (0) _ _ r (0) r (nf) _ _ r (nf) T (0) _ _ T (0) T (nf) _ _ T (nf) 

U P — 0|| , <-/p — 0|| , Jp j± — -t|| i± , -£p ± — -l|| i± , 



(f) 



-cf 



41n^ + 2L 



4 + 4^ 



-Cf 



4M B 
E(l-u)' 



T (f) 



0. (B.2) 



Here the shift of the 6-quark mass in O-j from the MS to the PS scheme has been taken 



into account in cjp . All quantities L,Cf, c{ 0) , cj nf) and 7] 



[9]. The expressions for Cp- 1 and Tp 1 ^ in (B.2) agree with — C|| tJ and — Tjj- 1 ^ given in [10], 
respectively, where a different definition of the longitudinal form factor £y with respect to 
the one used in previous works [9, 12] is employed. 



^°±,Tu n ± in (B.2) are given in 



(f) 



,(f) 
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